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The experimental data of BELLE Collaboration on the exclusive charmed meson 
pair production in the process of monophotonic e + e~-annihilation (e + e~ — > 7* — ► 
DD) has been studied. It has been shown that these data is described satisfactorily 
in the frame work of constituent quark model. Our studies have demonstrated that 
the central production process e + e~ — * e + e~77 — > e + e~ DD + X and the process of 
monophotonic e + e~-annihilation yield comparable numbers of the charmed meson 
pairs. 



The exclusive meson pair production is the unique opportunity to research the asymptotic 
behaviour of the meson form factor in the frame work of pQCD. The heavy meson pair production 
is of special interest. In this case the effective heavy quark theory allow us to connect the meson 
production processes with the meson decay processes [1]. The asymptotic form factor behaviour 
is described by the factorized amplitude which can be represented as a wave function f(x,Q) 
multiplied by a hard interaction amplitude [2]. In the leading order of a s only one gluon exchange 
contribute into the hard interaction amplitude and the wave function f(x,Q) have peak at x — 
A-qcd/M, where Aq C d is a strong interaction scale, and M is heavy quark mass. The distribution 
of f(x,Q) over x becomes thiner with increasing heavy quark mass. Due to this feature one can 
use the following approximation for the wave function: 



In this approach the momentum fraction carried out by quark is proportional to the quark 
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1. INTRODUCTION 
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mass: 
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The models of such kind [2] are taking into account the hard "tail" of the wave function 
only and neglect the soft component of the wave function. Nevertheless these models correctly 
reproduce the cross section behaviour near the threshould as well as asymptotic behaviour at 
Q 2 — > oo. Moreover, difference between the peak approximation calculation and the calculation 
beyond the peak approximation (1) is not essential [3]. In this work we endeavour to study 
the recent experimental results of BELLE collaboration [4] for D* + D*~ -pair production in e + e~- 
annihilation at y/s = 10.6 GeV. Also we predict the cross section value of the charmed meson pair 
production in 77-interaction at BELLE. 



2. £ + £--ANNIHILATION 

The following cross section values of the charmed meson pair production in e + e~ -annihilation 
have been measured by BELLE Collaboration^]: 

a(e + e- -> D* + D*-) = 0.65 ± 0.04 ± 0.07 pb, 

a(e + e- -> D + D*~) = 0.71 ± 0.05 ± 0.09 pb. 

The differential cross section distributions over azimuthal angle have been measured too. 

The interaction energy at BELLE (10.6 GeV) is large enough to determinate the leading asymp- 
totic contribution into the meson form factor: 
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where fu is the constant of leptonic decay, and Q 2 = s = (k e + + k e -) 2 is the e + e~-interaction 
energy squared. Nevertheless, it is worth to mention that the virtuality of gluon splited into 
quark-antiquark pair qq is slightly smaller than this energy 

Q ~ x q Q , 

where x q is the momentum fraction of the meson carried out by the light quark (See Fig. 1). If 
the pre-asymptotic terms would be large, then the model prediction would not reliable. But as it 
was shown in paper [3], the pre-asymptotic terms are not essential. 

Let us return to the discussion of the constituent quark model. As it mentioned above, the cross 
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section behaviour near the threshould is described correctly by this model: a ~ (| — M 2 ) , 
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where L is orbital momentum of the meson pair. (For production in e + e~-annihilation L = 1.) 
At M — > oo the constituent quark model, as well as the effective theory of heavy quarks give the 
same predictions for the cross section ratio: 

ffpp : crpv '■ Oyy = 1 : 4 : 7, 

where P and V are pseudoscalar state and vector state, correspondingly [5]. At large energy the 
cross sections <Jpp, opy and oyy behave as follows: 

a 2 a 2 s 
cr ~ — — . 

An angular distribution of the cross sections have the following simple form: 

da(DD) 



dcos 



sin' 6, 



da(DD*) 1 „ 
, V ' ~ 1 + cos 2 9, 
a cos 6 

where is the angle between the final meson direction and the initial beam direction in the 
center-of-mass system. 

For the case of DD* pair production only transversal component of D* contributes into cross 
section, because the production of longitudinal D* component do not allowed by lows of parity 
and angular momentum conservation. Indeed, in the center-of-mass system the matrix element of 
DD* production in 1~ state looks like follows: 

[e D .(p) x p] • e qV ?D, 

where e^. and e q are polarizations of D* and virtual photon, correspondingly. One can see from 
this formula, that the longitudinal component contribution D* 

en = n(ne D »j — , 

m 

equals to zero ( n = p/|p|, E and M are energy and mass of D* meson). 

The longitudinally polarized D* meson could be produced in the process of diphotonic e + e~- 
annihilation, but as is was shown in [6], the charmed meson yield is small for this process. 

The pair production of vector mesons D* and D* in the center-of-mass system can be describe 
by two independent structures: 



Mi ~ (eip)(e 2 e q ) + (e 2 p)(eie q ), 
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M 2 ~ (eie 2 )(pe q ). 

Both amplitudes correspond to 1~ state of the dimeson system. First amplitude contribute to 
D* L D?p-p&ir production. Second one contribute to all possible variants of particle polarizations. 
The ratio between these amplitudes depends on the quark masses. 

The satisfactory description of the experimental data of BELLE Collaboration has been achieved 
for the following values of the model parameters: 

a s = 0.3, 

f D = 200 MeV, 
m q = 0.17 GeV, 
m c = 1.5 GeV. 
We have predicted the following cross section values: 

a(D* + D*-) : a(D + D*~) : a(D + D~) = 0.73 pb : 0.58 pb : 0.02 pb 

One can see that the production of two pseudoscalar mesons is suppressed by the order of 
magnitude. 

The experimental angular cross section distribution da/dcosQ for the processes of D*D*- and 
DD*-paiv production in comparison with the model predictions has been shown in Fig. 2. The 
best description of the experimental data is achieved for m q = 0.17 GeV and m c = 1.5 GeV. 

It is worth to mention that the model under discussion is rather rough approach, which takes 
into account only hard "tail" of the meson wave function. The contribution of the soft part of 
wave function has been neglected. That is why the cross section dependence on energy have not 
minimum for the case of production of the two pseudocalar mesons, which has been predicted in 
the frame work of the effective theory of heavy quarks [1]. However cross section values predicted 
in [1] are large than experimental data by the order of magnitude. It is worth to mention ones 
more, that the results of calculations above the peak approximation (1) do not differ essentially 
from the results obtained in the frame work of peak approximation [3]. 

3. PHOTONIC PRODUCTION OF £)WdW-PAIR 

The exclusive production of D^D^-pair in the photon-photon interaction has been studied 
in our paper [7]. Our results have been compared with the predictions of the effective theory of 
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heavy quarks. In the constituent quark approach the photonic production of the charmed meson 
pair is described by twenty Feynman tree level diagrams (see Fig. 3). The analogous diagram set 
has been used to calculate the inclusive P c -meson production cross section [8]. This diagram set 
can be subdivided into three gauge invariant parts. The first group of diagrams (1-6 of Fig. 3) 
corresponds to the case when heavy quark radiates gluon which splits into light quark pair. The 
second group ( diagrams 7-12) can be received from the first one by permutation Q <-> q. The 
third diagram group (diagrams 13-20 of Fig. 3) corresponds to the independent production of the 
of the quark pairs ( 7 — > QQ, 7 — > qq) followed by there fusion into mesons. 

The third diagram group dominates at relatively small transverse momenta. Due to the con- 
tribution of such diagrams the factorization theorem can not be applied to describe the inclusive 
production of P c -meson in the wide kinematic region [8]. In particular, the strong cross section de- 
pendence on light quark charge e q exists near the threshould as well as at large interaction energy. 
Also, we have observed such e q dependence for the process of photonic D^D^ pair production. 
That process feature has underlined also in the paper [9] for the region of large energies. 

In this article we continue our study of the constituent quark model prediction for different 
kinematic regions. Our analysis is based on diagrams of Fig. 3. 

Let us consider the I? (*)/}(*) pair production near the threshould. As it has been mentioned 
above, D*D* pairs and DD ones are produced near the threshould in ^-wave state. (The cross 
section depends on momentum k between P-mesons in the center-of-mass system as ~ k 2L+1 ). 
This state corresponds to the total momentum and parity of the diphoton system + . On the 
contrary, P*P-pair is produced in P-wave, because S'-wave state production is not allowed by the 
Landau- Yang theorem. 

5- wave state of the meson pair is produced from the 77 system 1 + . P-wave state of the meson 
pair D*D corresponds to the 77 system 0~, thus the cross section of the D*D production increases 
near threshold as k 3 . As one can see in Fig. 4, the cross section values of P^P^-pair production 
depend on light quark charge, whereas quantum numbers of the meson pair, as well as energy 
behaviour of the cross sections do not depend on this charge value. Also, one can see, that the 
cross section distribution over angle for neutral meson pair and charged meson pair differ from 
each other considerably. 

Our model predicts practically isotropic production of the P°P -, P*°P -, P°P* -, D + D~- and 
P* + P*~-pairs near the threshould: 

da 

— const, 

a cos 6 
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whereas D + D*~ -pair produce peripherally even near the threshould (see Fig. 5). However, as it 
can be clearly seen in Fig. 4, the production of such pairs is suppressed near the threshould. It is 
worth to mention, that the total cross section value of the neutral meson production is three times 
large than the cross section value of the charged meson production. 

At large energies the cross section distributions over cosG have peripheral form with striking 
maxima at +1 and -1. This fact has been also emphasized in the paper [9], where the analytic 
form of amplitude of the pseudocalar meson pair production has been performed at hight energy 
limit: 



2 1 + p 2 cosO 2 2 

6q} i - p cos e 2 + 6q 



where 



s 

At large value of s one can see the difference between the asymptotic behaviour of the cross sections 
for the charged and neutral meson production. The production cross section depends on s as 1/s 2 
for all charged meson pair. On the contrary, the production cross sections for neutral meson pairs 
have different asymptotes: 

1 1 1 

<T D o D o : <J D *o D o : <J D *o D *o ~ — : — : — . 

S S S 

However, it is worth to mention that the cross section behaviour becomes asymptotic after 
y/s > 20 GeV, therefore the study of it have the theoretical importance only. In addition, the next 
order logarithmic corrections can play essential role at such energies. 

4. £)W^W-PAIR PRODUCTION IN 77-INTERACTION AT BELLE 
It has been shown in the paper [6] that the contribution of diphotonic annihilation 

e+e" -> 7*7* -> DD (2) 

into the total cross section is small. 

Nevertheless, another opportunity to produce charmed meson pairs in the photon-photon in- 
teraction exists. Charmed mesons could be produced in the interaction between the effective 
(equivalent) photons radiated by the initial fermions: 

e + e~ — > e + 7e _ 7 — > DD + e + + e~. (3) 
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The cross section of this process is described in the frame work of the "partonic" model: 

cr = J J a^(s)f / (x 1 )P(x 2 )dx 1 dx 2 , 

where X\ and x 2 are momentum fractions of the initial electron and positron, correspondingly, 
carried out by interacting photons. The photonic density P{x) obeys Weizecker- Williams formula: 

p(x) = f(\ - (1 xf) In ^™ ax m\x 2 ( — —Y) » 

27TX y Qmin V Qmin Qmax / / 

where m e is the mass of electron, Q min = m 2 e x 2 /(l — x), and Q max — 1 GeV 2 . 

In contrast to the diphotonic annihilation (2) the amplitude of process (3) is not suppressed 
by s-canal propagators. Thus it would be interesting to estimate the cross section value for the 
process of the charmed meson production in the interaction between the effective photons. 

For the parameter values used to estimate the cross section in e + e~-annihilation we have 
received the following unexpected results for the process (3): 

a(D* + D*-) : a(D + D*~) : a(D + D~) = 1.52 pb : 0.33 pb : 0.13 pb. 

a(D*°D*°) : a(D°D*°) : a(D°D°) = 1.39 pb : 0.43 pb : 0.40 pb. 

It is clear that these cross section values are comparable with ones for the process of e + e~- 
annihilation. One can see in Fig. 6, that the distribution shapes do not differ essentially from ones 
for the case of e + e~-annihilation. It is obvious, that the meson energy is less than the energy of 
meson produced in e + e~-annihilation. The distribution over z = ^=^- is shown in Fig. 7, where 
Pd is the D meson momentum in the center-of-mass system of e + e~-pair. One can see from that 
distribution that the averaged value of z is about 0.3 for the such process. Therefore the averaged 
meson momentum is about 1 GeV. 

It is worth to notice that the model under discussion predicts the dominance of the vector- vector 
pair production. 

CONCLUSION 

We have shown that the experimental data on charmed meson pair production at BELLE 
(a/s = 10.6 GeV) is described satisfactorily in the frame work of the constituent quark model. In 
that approach the hard part of the amplitude can be calculated with the help of pQCD motivated 
diagrams. 
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The production mechanism in e + e~-annihilation is rather simple. The virtual photon splits into 
cc-pair followed by hardonization process. Thus the cross section value do not depend on light 
quark charge. 

In contrast to e + e~-annihilation in 77-interaction the cross section values depends strongly on 
the light quark charge at low energy, as well as at high energy. The essence of matter is that in 
the process 77 — > DD the light quark interaction with the photonic field can not be neglected. 

It has been shown that the cross section value of the central D meson production in the process 

e + e~ — > e + e _ 77 — > e + e~ DD + X 
and the cross section value of D meson production in e + e~ -annihilation are comparable. 
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Figure 1. Feynman diagrams for the process of charmed meson pair production in e e + -annihilation. 
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cos 

Figure 2. The cross section distribution over cos 8 (pb) in comparison with the experimental data for the 
process e+e~ -> D* D* (a) and for the process e+e" -> (b). m c = 1.5 GeV, m, = 0.17 GeV, f D = 200 

MeV, a s = 0.3. 



1 2 3 




10 11 12 




13 U 15 16 




17 18 



Figure 3. Feynman diagrams for the charmed meson pair production in the photon-photon interaction. 
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Figure 4. The cross sections app (o, •), o PV (A, ▲), oyy (□, ■) as a function of k = ^Js — s t h for neutral 
Z?-mesons (closed symbols) and for charged £>-mesons (opened symbols). 
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Figure 5. The cross section distribution over cos for the processes of D + D* and D D* + pair production 



in the H e + e -annihilation at k = y/s — s t h = 0.1 GeV. (Other meson pairs are produced at small k 
isotropically.) 
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Figure 6. The differential cross section dayv/d cos 9 (solid curve), dupyjd cos (dashed curve) and 
dapp/dcosQ (dotted curve) for the processes of charged meson production (a) and neutral meson production 
in the effective photons interaction. 
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Figure 7. The cross section distributions over z = for the processes of charged meson production (a) 

and neutral meson production (b) in the effective photon interaction. The designations are the same as in 
Fig. 6. 



